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Abstract 

The rate of step propagation for an isolated, continuous and curved step with bulk diffusion is solved within 
the steady-state approximation and it is noted that the Chernov model of parallel sequence of step advance 
is without effort adapted to the curved step. The solution of the Laplace equation in toroidal coordinates 
for a diffusion layer of arbitrary extension yields the diffusion flux density to the step perimeter expressed 
as a Legendre series. For the special case of diffusion without advection (infinitely thick diffusion layer), the 
step velocity is found to vanish as an inverse square-root of the radius of curvature for large values of the 
same. For a sequence of steps, this effect should be even more pronounced. 
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1. Introduction 

While it can be argued that between surface (2D) and bulk (3D) diffusion there is no clear-cut physical 
distinction, the former being geometrically contained in the latter, the division has found great practical 
value in the mathematical treatment. Burton, Cabrera and Frank [1[ solved the surface-diffusion equations 
for mass transport to steps. Chernov [2| was the first to solve the case of pure bulk diffusion to the straight 
step, in isolation or in parallel sequence. This was later generalized by Gilmer, Chez and Cabrera Q to 
account also for a diffusion of growth units to the surface from bulk. The coupled volume-surface diffusion 
problem was then attacked admirably by van der Eerden [4[ using a method where coupled integral equations 
define the total flux density due both to surface and bulk contributions. Van der Eerden obtained explicit 
solutions for a train of parallel, straight steps and for an isolated curved step. The unacquainted reader is 
referred to these previous works for the necessary background. 

As pertains to growth from solution, these treatments all assume a concentration gradient of finite 
extension separating the bulk concentration from that at the steps: the so-dubbed "unstirred boundary 
layer." While in cases of growth from stirred solution, this is reasonable, in cases where advection is kept 
to a minimum, this assumption can be dispensed of if the curvature of the steps is properly taken into 
account. Especially as in some cases of growth from solution, stirring may be undesirable [5]. In this Paper 
we solve the bulk diffusion problem up to an isolated curved step for an arbitrarily thick diffusion layer and 
also extend Chernov's treatment [3] of parallel steps to the case of an infinite sequence of parallel, curved 
steps of equal curvature. Equations are obtained where the rate of step advance is expressed as individual 
contributions from each angle of impingement at the step, so that the total step velocity is given as an 
integral over the surface of a semitorus. It will be shown that, the Gibbs-Thomson effect notwithstanding, 
under bulk diffusion from infinitely far away the step velocity should decrease after an early maximum. The 
main assumptions are 



• The crystal surface is mathematically flat. 
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• The step is a segment of a semitorus of greater radius R and smaller radius a/V, where a is the 
molecular diameter. 

• The step velocity is very slow compared to that of the incoming growth units. 

• Activity coefficients are unity. 

• There is a steady-state. 

• Surface diffusion is negligible. 

• There is no Ehrlich-Schwoebel barrier. 

• Kinks are so numerous along the step that it is essentially a continuous sink for growth units. 
Of these, the last three may conceivably be relaxed, at the expense of added parameters. 

2. Toroidal coordinates and the solution of the Laplace equation 

The symmetry of a curved step suggests the use of toroidal coordinates. These are obtained by rotating 
the 2D bipolar coordinates around a fixed axis of symmetry and may be defined as @, p. 112] 

r sinh ri cos ib 

x = u a 1 

cosh n — cos 9 

r sinh rj sin ip 
^ cosh rj — cos 9 

z = r J^l (3) 

cosh n — cos 9 

where the coordinates satisfy 77 £ [0, 00), 9 G [— tt,tt] and ip € [— 7r, it], respectively. The parameter r is a 
scaling factor of the coordinate system. For each value of r, a new toroidal coordinate system is obtained. 
It is related to the greater radius of the torus, R, by 

r = i?tanhryi (4) 

and to the lesser radius a/ir by 

r = — sinh 7/1 (5) 

where 771 defines the surface of the torus. The angle 9 refers to the latitudinal direction, where 9 = ir/2 is 
directed parallel to the z-axis. The azimuthal angle ip denotes rotation about the z-axis at a distance 

r sinh rj 



cosh rj — cos 9 



2.1. General solution of the Laplace equation 

The Laplacian in toroidal coordinates is given by @, p. 112] 



y2= (£-cos0) 3 



£-cos# / \4 — cosf J £ — cos 9 



(6) 



where £ = cosh 77. Denoting the concentration function by c(^,9,il)), the Laplace equation is 

V 2 c(£,0,V)=O (7) 



This equation is separable and if c(£, 8, ip) is expressed as [6j, pp. 112-115] 



c(£, 6,ijj) = A + — cos 6H{>q)®{e)*{il)) (8) 
where A is a constant, the solutions are obtained from 

(£ 2 - l)H" + 2£ff' - {{p 2 - 1/4) + q 2 /{e - l))H = (9) 

Q"+ P 2 e = o (io) 

*" + 9 2 * = (11) 

where the primes indicate derivation. The first of these is the associated Legendre equation (p, q are con- 
stants) and its solutions are linear combinations of the associated Legendre functions [3] {P q _i(0} an d 

P 2 

{Q q _i {(,)}■ Of these, the Q-functions become irregular for r) — ► 0, corresponding to great distance from the 

P 2 

torus or close to the z-axis and must be excluded for physical reasons. The solutions to Eqs (flT))) and pip 
are of the same general form. Because of the requirement that the solutions be periodic on the torus, there 
is no linear term even when p = q = 0, so they are 

0(0) = B p cos( P 9) + C p sm(p9) (12) 
*(V) = B' q cos(#) + C' q sin(#) (13) 

andp, q are non-negative integers because of the same periodicity requirement, whereas B p , B' q and C p , C' q are 
arbitrary constants. Assuming azimuthal symmetry, Legendre's equation replaces the associated Legendre 
equation, B' q = C' q = and the general solution is given by 

oo 

c(e, 9) = A + cos (B p cos(p0) + C p sin(p0))P J) _i(O (14) 

p=0 

36 Particular solution for a continuous step at arbitrary distance from bulk in radial direction 

The coordinate £1 defines the surface of the step where the concentration is 

c(e 1 ,0)=c 1 (0) (15) 

whereas £5 < £1 defines a toroidal surface at a distance 6 from the step (the unstirred boundary layer), 
where the concentration has reached its bulk value, 

c(t5,6)= Coo (16) 

This means that 

00 

Cl (0) =A+y/£i- cos J2( B p C0S (P 9 ) + C p sm(p0))P p _ 1 (&) (17) 

and 

00 

coo = A + - cos0^(B p cos(p0) + C p sin(p0))P p _i(^) (18) 

Two boundary conditions and three classes of unknowns (A, B p , and C p ), we need a further condition and 
impose the physically reasonable case that Coo is even in 8, in which case all sine terms must vanish. Then, 
through the relation (Heine's equation) 

1 /2 00 

a = — W P -|(£)cos(p0) (19) 

VC - COS 7T ^ 2 
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we may deduce that 

V2( Cl -A)(2-So p ) Q P -j(£i) 
b p = ^ ~p fTT ( 20 J 



and 



N /2(c oo -A)(2-5 0p ) ( 9j>-ife) 



7T 



Summing both equations and then solving for A, yields 

siKi)-Ei(e«) 



where we have defined _ 

7T P p -l(x 



p=0 

for future reference. Finally, 

£i = i?7r/a (24) 
and, if 5 << i?, i. e. for a sufficiently great radius of curvature, 

Cs = R/S (25) 

otherwise, Vx, y, z\x 2 + y 2 + z 2 = S 2 , we have 



m (^Tg±l)j±i! (26) 

V \Wx 2 +y 2 - r) 2 + z 2 ) ) 

by virtue of the inverse coordinate transformation. It is to be noted from this last expression that £s — > 1 
when (5 — > oo. 

Equidistant parallel steps: Chernov solution 
In Chernov's solution for the straight step the solution for an isolated step is generalized to the 
sequence of parallel steps by substituting the radial distance by a periodic term. On the condition that all 
the steps have the same radius of curvature, are weakly curved and/or at large separation, this approach 
can be applied also the current case. Making use of the expression for £ valid when R is great, Chernov's 
change of variables takes the following form, 

^ 1 = (27) 

/sin^Wsinh 2 (^) 



Here the steps are spaced d apart, and — \/ x 2 + y 2 is the coordinate which runs over their separation, 
whereas a — z is the surface normal coordinate. The boundary conditions are now (see |2j|) 

lim c = a (28) 

a — >a/ 7r 

lim c = Coo (29) 

3—0 
a — >6 

The boundary conditions affect the following changes 

( ° - 4?) (31) 

but the equations are otherwise unchanged. For future convenience, we also define p = a 2 + 1 . 
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3. Steady-state step velocity 

In toroidal coordinates, the gradient operator is given by 



cosh 77 — cos 8 Ud if) d\ 

r I drj dO sinh rj dip J 



m = _D(™^-™°) (33) 



so the diffusion flux density to the embryo perimeter, as given by Fick's law, is 

'dc 
dij 

In the general case, the step velocity is given by § 

Voo{R) = v m {a 2 /x )v{ci - co)(l - R*/R) (34) 

42 where v m is the molecular volume of the building unit in the crystal, v an attachment rate constant, 

43 ci the concentration of growth units in the step vicinity, cq the solubility, R the radius of the surface 

44 cluster, R* the critical radius for surface nucleation and xq the mean distance between kinks. The factor 

45 (a 2 / xq)v(c\ — Co) = ji is the integration flux density. Balancing the integration and diffusion flux densities, 

46 means that the former becomes implicitly dependent on 8 and that the step velocity can be decomposed 

47 into contributions from flux densities at different angles. We shall treat the different cases individually. 

48 3.1. Step at arbitrary distance from bulk in radial direction 

In this case, 

m = (^^E 2 + & - ooefl)*^) (35) 

since = sinh 771 = irr/a and where 

00 

Z 2 =J2B P cos(p6)P p _i(Z 1 ) (36) 

p=0 

and = d'E^/d^x. Not surprising, the flux density is higher from inside the perimeter, as we have neglected 
the Ehrlich-Schwoebel effect [9|, [lOj • Solving ji = jd for c\ (8) and inserting it into Eq. ([3"4"j) gives 



a 2 Ei(^)(coo - co)Dnv m vy/Ci-cos0(H 3 + 2E^i - 2E 3 cos8) (l - 
V °° = 2a 3 (E 1 (^) - Ei + £1(^)^0 v / ei^o¥0(E 3 + 2E 3 a - 2E 3 cos0) 

Here E 3 = (y/2/%) £~ (2 ~ M cos(pfl)P p _ 1 (6) and E£ = 0E 3 /0fi. 



(37) 



50 3.2. Chernov-type solution 

In this case, it is rather jd = — Z?(dc/p) p= » that we seek. We do not have an explicit expression for c, 
but we may write by virtue of the chain rule 



dc £ — cos 8 dc ( drj da drj dj3 



dp r dn \da dp dj3 dp 



(38) 



where we have added a correction factor 4> whose purpose is to make the overlapped equation coincide with 
the isolated one, in the limit of infinite step separation. In this expression, the first factor is recognized as the 
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non-overlap adjusted gradient, so the bracketed factor should tend to 1/ 4> when d — > oo. Now da/dp = p/a 
and dp /dp = p/p from the definition of p, whereas 

dri _ 7ri?cosh(^)sinh(^) 
da 



d(sin 2 (^) + smh 2 (^) 



(39) 



R? 



Ksinh ( 



- 1 



and 



drj 
~dp~ 



ttR cos 



(f)-(f) 



rf(^ 2 (f)+sinh 2 (^ 



/,'- 



- 1 



from the substitution of variables, Eq. (f2"T|) . In all, 

/ 



M(0) = 



ttD ( a/^i — cos ( 



S 2 + -cos0)*E£ 



(40) 



\ 



V d ( sin2 W + sinh2 »)) 3 - v 

(41) 

We shall denote the last factor by —f(0) for brevity and where the minus sign is for later aesthetic reasons. 
Its dependence on 8 stems from all possible a and /3, such that a 2 + /3 2 = (a/ir) 2 . They are given by 
a\ = (a/7r)sin(9 and /3i = (a/n) cos6. Solving the flux density balance for c\{9) and inserting it into Woo 
yields 



a 2 Si(6)(coo - co)£>/(0)7™ m iV$i-cos0(£ 3 + 2E' 3 & - 2E^ cos0) (l - ^) 
2a3(E 1 fe) - Ei(&))i/ + £i(ft)£>/(0)7ra; o V£i - cos0(E 3 + 2E 3 £i 1 2£ 3 cos0) 



(42) 



4. Non-advective step velocity 

While it is not expected that Eq. (f37|) capture all the fine peculiarities of the rate dependence on stirring 
of growth from solutions, we do note that it is well-defined for all possible extensions of the stagnant layer. 
This parameter plays a role in the balance between diffusion and kinetic control of the growth. One case that 
is of particular interest, since it has to the author's knowledge not been treated before, is that of complete 
absence of advection. In this case we let 8 — > oo and we have 



a 2 (coo - c )£>7ro m iV£i - cos (9 (E 3 + 2£ 3 £i - 2E 3 cos(9) (l - 



(43) 



2a 3 v + Dirxo V£i - cos 9 (E 3 + 2E 3 £i - 2E 3 cos( 
Because of the asymptotic properties of the Legendre functions, Q 

E 2 ex 1 (44) 

for large R, so for this expression, when R — > oo, the step velocity vanishes as 0(R~i). In the case 
of growth due to surface nucleation, [ll|, [Til, [l3| especially in the mononuclear regime, this means that 
because the surface covered (from a surface cluster) grows as R 2 , however, the surface covered actually goes 
asymptotically as 0(t^) (t is time). On the other hand, in the case of complete advection, £5 = £1 and the 
expression reduces to 

_ a 2 (coo - co)v m v ( R* 



\ 1 - IT (45) 

x \ R J 

which attains its maximum when R — > 00 and the surface coverage goes like 0(t 2 ). As pertains to Eq. (|42p. 
compared to Eq. Q37p . the overlap correction is completely contained in —/(#), which is a monotonically 
decreasing function of R. In spiral growth [l|, where such an equation would find its use, the effect 
should hence be even more pronounced. This mass transport restriction on the rate of step advance works 
to counter the Gibbs-Thomson increase. 
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5. Conclusion 



Within the steady-state approximation, a general equation for the isolated step velocity with bulk dif- 
fusion from arbitrary distance, taking into account also the inhomogeneity of the resulting diffusion-field 
around it, has been derived and the method of Chernov Q has been adjusted to take account of that correc- 
tion so that a sequence of equally curved steps can be described. For the special case of bulk diffusion from 
infinite distance (no advection) , it has been shown that the step velocity for great radii of curvature should 
vary as an inverse square-root of the same and that this effect should be less pronounced for an isolated step 
front, than for a sequence of many. 
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